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Abstract 

We use Monte Carlo methods to directly evaluate D-dimensional SU(N) 
Yang- Mills partition functions reduced to zero Euclidean dimensions, with 
and without super symmetry. In the non-supersymmetric case, we find that 
the integrals exist for D = 3, ./V > 4 and D = 4, N > 3 and, lastly, D > 5, 
N > 2. We conclude that the D = 3 and D = 4 integrals exist in the large N 
limit, and therefore lead to a well-defined, new type of Eguchi-Kawai reduced 
gauge theory. For the supersymmetric case, we check, up to 377(5), recently 
proposed exact formulas for the D = 4 and D = 6 D-instanton integrals, 
including the explicit form of the normalization factor needed to interpret the 
integrals as the bulk contribution to the Witten index. 
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In a recent study we developed Monte Carlo methods in order to deal with di- 
mensionally reduced Yang-Mills theories. We directly work with the gauge potentials, in 
contradistinction to more conventional numerical studies of lattice gauge theories. Our ini- 
tial interest was in establishing reliable methods for the numerical calculation of the bulk 
contribution to the Witten index in supersymmetric field theories. Some further results in 
this direction are presented below. The methods are actually even more efficient if applied to 
the non-supersymmetric case, and allow us to settle the question of existence of Yang-Mills 
theory dimensionally reduced to zero dimensions in a surprising way. 

Let us first consider maximally reduced D-dimensional SU(N) Euclidean Yang-Mills 
theory. The formal functional integral for the partition function then becomes an ordinary 
integral: 
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Since there are directions in the integration space which are not suppressed, these integrals 
were generally believed to be ill-defined. Indeed, it is very easy to show that e.g. for D = 2 
the integral diverges for all N. Nevertheless, in we obtained an astonishing result: For 
gauge group SU{2) the exact result reads 
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Therefore, the reduced bosonic theory is not necessarily divergent. It is natural to ask how 
eq.flU) generalizes if iV > 3. The methods used to derive eq.(|]) are specific to SU(2), and 
no analytic result is known for higher gauge groups with N > 3, except for D = 2, as 
mentioned. However, we can modify the methods of Jl| to decide the question of existence 
by Monte Carlo evaluation. Our results suggest the following intriguing answer: 
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In particular, this suggests that a well defined large iV limit exists in dimensions D = 3 
and, most interestingly, D — 4. It opens the exciting possibility that appropriate large iV 
correlation functions computed for the model (|I|) can be related to large iV Yang-Mills field 
theory through the Eguchi-Kawai mechanism [Q. E.g. , one could consider Wilson loop 
operators such as 
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in the limit N — > oo. In fact, our model eq.([l|) is a new type of continuum (as opposed 
to lattice) Eguchi-Kawai model. Reduced models are frequently plagued by the need to 
introduce quenching, but our models are already in the weak coupling phase, so this does 
not appear to be a problem. However, these important questions are beyond the scope of 
the present analysis. 

Let us explain how to obtain convincing evidence for the result eq.(§). In U we computed 
ratios of integrals for different dimensions D. Such a strategy is rendered necessary by the 
strong fluctuations of the integrand in eq.(|l|). One key point in our procedure [[| consists 
in compactifying the integrals: after introducing polar coordinates (xi,...,Xd) = (£ld,R) 
(d = D(N 2 — 1) being the total dimension of the integral), we exactly perform the R- 
integration. In the present (bosonic) case of eq(P, Zd,n can be written as 
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To obtain the absolute value of Z^ jy, we now consider a series of interpolating functions 



(fid) with 1 = a < Oil < 



< ai-i < ai = (notice that zd,n > 0). These 



interpolating functions allow us to connect zd,n{Q<i) to a constant in much the same way as 
we can (sometimes) compute the free energy of a statistical physics system by simulating at 
various temperatures from oo down to the temperature of interest. The term with ai = 
plays the role of the exactly solvable high-temperature limit, since we can integrate the 

1 analytically on the <i-sphere. We now write 
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Each of the terms [ ] in eq(|7|) is computed in a separate Monte Carlo run, in which points are 
sampled according to 7r(0) ~ z a \ using the Metropolis algorithm. The fluctuations of the 
operator to be evaluated, z ai -^~ a \ are much damped if ~ a^. In practice, we have used 
I = 10, corresponding to the number of work stations available to us. Judiciously chosen 
values of ctj, closely spaced as we approach 1, allow us to compute any finite bosonic integral 
with ease. As a simple check, we have computed Z n=2j d=5 ( c / ec L-(@)) to within 1% precision 
in about 30 minutes of individual CPU time on 10 machines. 

The finiteness of the integrals was also checked with the extremely powerful qualitative 
Monte Carlo method described in Q . In that method, we solely monitor the autocorrelation 
function of the Metropolis random walk. The cases analytically known to diverge were 
easily identified (D = 2 for all N, D = 2, 3, 4 for N = 2). Of these cases, the N = 2, D = A 
random walk appears less divergent, which agrees with the analytic result that the divergence 
is marginal (the N = 2,D = A + e integral exists). A similar behavior is observed for 
iV = 3, D = 3, which we believe to be marginally divergent as well. 

Using the techniques just described, we are able to present the following table: 
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Table V. Direct evaluation of bosonic Yang-Mills integrals 
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The computation becomes simpler both with D (it is completely trivial for D > 5) and 
N. We are thus confident about the statements expressed for large N and arrive at the 
prediction @ . This is the central observation of the present paper. 

Turning now to the supersymmetric case, the integrals reduced to zero dimensions read 
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The only possible dimensions are D = 3,4, 6, 10 corresponding to M = 2, 4, 8, 16 real super- 
symmetries, respectively. Integrating out the fermions, we find an integral that differs from 
eq.(|ID only by a modified measure: 
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The integrand is weighted by a very special homogeneous polynomial Vd,n of degree k = 
(D - 2)(N 2 - 1) = \M(N 2 - 1) in the variables X*. See e.g. [§ for further details. In a 
very recent calculation, Moore et.al. 0, heavily using cohomological field theory techniques 
(and therefore supersymmetry), evaluated, for D = 4,6,10 analytically a class of related 
integrals which are initially defined in Minkowski space and are Wick rotated in the course 
of the evaluation. These cohomologically deformed integrals are argued, after an appropriate 
prescription for the Wick rotation is found, to take the same values as the integrals (||). They 
find the integrals to be of the form 
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Let us point out that the D — 10 result had previously been conjectured by Green and 
Gutperle 0. For SU(3), we had also performed a careful numerical check of the form of 
eq.flTOD in J 7 ^ is a group-theoretic factor not worked out explicitely in 0]. In fact, it is 
unclear to us how to derive the the explicit form of J 7 ^ from their approach. In Q, it was 
found to be 
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It is important to stress again our finding that the Euclidean super symmetric partition 
sums are not some formal mathematical entities, but perfectly converging ordinary 

multiple integrals, just as their bosonic counterparts. It is very tempting to speculate here 
once more that appropriate correlation functions computed for N — > oo could be related to 
correlators of the full susy gauge field theories. 

For completeness, let us sketch the derivation of (|H|), in particular since its precise 
relationship with the non-explicit normalization of remains obscure. In || is was argued 
how the D-dimensional Euclidean matrix model emerges from the functional integral for the 
super symmetric gauge quantum mechanics of D — 1 matrices when computing the bulk part 
of the Witten index 

lim /3 _ Tr(-l) F e- /3H (12) 

One needs to integrate over the group SU(N) in order to project onto gauge invariant states. 
This integration becomes, in the limit (3 — > 0, an integration over the hermitean generators 
of the group 

1 JV2 ~ 1 rlY A 

vu~^n% us) 

•*N a=1 V^TT 

where Tn is precisely the constant relating the bulk part of the index to the Euclidean 
matrix model. It is easily found if we investigate how the normalized Haar measure T>U 
and the flat measure IIa^ 1 dX^ act on class functions. The latter depend only on the 
eigenvalues = e tXi of the unimodular unitary matrices U. One easily checks (e.g. by 
verifying orthonormality of group characters) that group integration on class functions f(U) 
is performed as 

J VUf(U) = ^ j fl ^7 A (^) A ( <)27T^(A 1 + • • • X N )f( Zl , ...,z N ) (14) 
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where 5 p is the 27r-periodic ^-function and A(^) = Tii<j( z i — z j)- On the other hand, the 
flat measure for the hermitean matrices Xp becomes 



N 2 -l jyi N n 

II ~i= = TT^^ II ^ 2 (^S{X 1 + ... + X N ) (15) 
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where the factors are easily checked by computing a Gaussian integral normalized to one 
(the Aj are the eigenvalues of Xr>). Then, replacing (since (3 — > 0, see 0) 

Zi = l + iXi + ..., (16) 

comparing ( p!4]) and ( p!5[ ) and multiplying by an extra factor of N, due to the fact that the 
unitary matrices localize on iV values because of invariance under the center of the group, 
we can read off from eq. flilf ) the result eq. (|TTD for Tn- 
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In order to test the group factor flTT|) as well as the form (|1~0"D of the D = 4 bulk part of 
the Witten index, we have applied the methods of direct Monte Carlo evaluation explained 
above to the D = 4 and N < 5 integrals. As conjectured in [1[], we find that the integrand 
is always positive if D = 4. This means that we can use the method of direct evaluation 
described earlier. The results are presented in table 2. 

Table 2: Direct evaluation of the D = 4 D-instanton integral 



N 


eqs.©,© 


Monte Carlo 


error 


2 


1.25 


1.25 


<0.01 


3 


3.22 


3.22 


<0.01 


4 


7.42 


7.6 


± 0.2 


5 


10.04 


10.2 


± 0.2 



We have also checked that the D = 6 integrals equal the D = 4 integrals to good precision 
by the ratio method of [[I]]. 

Finally, let us mention the open problem of the evaluation of the supersymmetric D = 3 
integral. In [flj, we conjectured 

Z^In = (17) 

This result is trivial for N even, but non-trivial for odd N. We suspect that the supersym- 
metric D = 3 integrals are absolutely convergent (except for SU(2), and possibly SU(3)), 
just like their non-supersymmetric counterparts, c/eq.(||), and that their well-defined value 
is given by eq.([H]) ||. Another interesting conjecture was presented in [[|, where it was 
suggested that a modified D = 3 integral Zd=£n, where the action is complemented by a 
Chern-Simons term, leads to 

Zd=In = jp- (18) 

Actually, for the solvable case of SU(2), it can be shown that eq.(|18D is not true, since one 
easily finds Z^^ N=2 = oo. However, our bosonic result (§) suggests that the susy integrals 
modified by Chern-Simons exist for at least iV > 4, and therefore it is feasible to test the 
conjecture eq. (|l8|) for generic gauge groups by our approach ||. 
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